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Summary 

'  Combinatorial  methods  are  used  to  derive  the  distribution  of 
the  busy  period  in  a  Markovian  tandem  queue  and  queues  with  finite  ca¬ 
pacity  in  five  models  each  involving  batches. 

V 

1 .  Introduction 

Tak^cs  in  his  book  (see  Tak/cs(1967))  and  many  of  his  papers 
amply  demonstrated  his  pioneering  work  on  the  application  of  combina¬ 
torial  methods  for  deriving  the  distribution  of  the  busy  period.  He 
(see  Takacs  (1962))  gave  a  new  direction  of  approaching  the  problem 
through  the  so-called  urn  problem  (which  is  a  generalization  of  the 
ballot  problem) .  The  technique  basically  requires  the  following  two 
steps: 

(i)  Write  down  a  finite  set  of  relations  in  the  form  of  inequalities  and 
equalities  on  a  sequence  of  random  variables  which  a  busy  period  when  suit¬ 
ably  conditioned  must  satisfy.  The  conditions  are  such  that  each  sequence 
has  the  same  weight  (in  terms  of  either  probability  mass  or  probability  density) 

(ii)  Count  the  number  of  possibilities  or  find  the  measure  of  the  set 
satisfying  the  inequalities. 

The  random  variables  are  in  terms  of  the  number  of  arrivals  or 

departures  or  the  arrival  or  departure  instants.  The  Markovian  assumption 

of  either  interarrival  times  or  service  times  plays  the  central  role  in 

>  *The  author  is  on  sabbatical  leave  from  McMaster  University  and  spending 
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the  combinatorial  technique  in  the  sense  that  it  assures  the  requirement 
of  uniformity  as  suggested  in  (i)  and  thus  leads  to  the  enumeration 
problem  in  (ii) . 

The  purpose  of  this  paper  is  to  further  explore  the  scope  of 
applications  of  the  same  technique.  With  this  objective  in  mind,  in 
Section  2  we  present  three  enumerative  propositions  which  find 
applications  in  the  remaining  part  of  the  paper. 

Although  queueing  network  models  abound  in  applications  (such 
as  communication  networks,  computer  time  sharing  and  multiprogramming 
systems)  the  study  of  such  systems  has  not  progressed  at  a  pace  conraen- 
surate  with  its  importance  (see  Jackson  (1957) ,  Jackson  (1963) , 

Disney  (1980)).  In  Section  3,  we  consider  the  simplest  form  of  a  queueing 
network,  viz.,  a  Markovian  queue  in  series  (also  called  a  tandem  queue) 
and  obtain  the  joint  distribution  of  the  length  of  a  busy  period  and  the 
number  of  arrivals. 

Queues  with  finite  capacity  is  rather  common  in  practice.  How¬ 
ever,  the  literature  on  busy  period  with  finite  capacity  is  not  large 
(e.g.  Cohen  (1971),  Enns  (1969),  Neuts  (1964),  Mohanty  (1972), 

Takacs  (1976)).  The  direct  combinatorial  approach  worked  well  on  de¬ 
termining  the  distribution  of  a  busy  period  with  infinite  capacity. 

Never- the-less,  it  has  limited  success  for  queues  with  finite  capacity 
(Mohanty  (1972)).  Not  surprisingly,  most  often  the  finite  capacity  case 
is  treated  through  the  Laplace-Steiltjes  transform  (Cohen  (1971), 

Enns  (1969),  Takacs  (1976))  rather  than  the  exact  distribution  with  the 
help  of  the  combinatorial  method  because  part  (ii)  evaluation  in  general 
is  not  simple.  However,  the  method  can  be  applied  to  more  situations 
than  in  the  past  and  in  Section  4  and  Section  5  we  obtain  the  disbribution 


of  a  busy  period  with  finite  capacity  (i.e.,  fixed  maximum  queue  length 
or  maximum  work  load)  in  four  different  models,  each  involving  batches. 


2.  Some  Auxiliary  Combinatorial  Results 

In  .this  section  we  present  three  combinatorial  results  each 
of  which  is  applied  in  the  subsequent  sections. 

Proposition  1.  (Kreweras  (1965),  p.  35) 

Let  N(k;  ja,  b)  denote  the  number  of  k-dimensional  lattice 
paths  from  the  point  a  ■  (a^,..., 

such  that  every  lattice  point  (x^ . x^)  on  the  path  satisfies  the 


a^)  to  the  point  t>  *  (b^,...,b^) 


condition 


X1  t  •**  t  V 


Then 


(1)  N(k;  a,  b)  -  ft  -  a] !  |  |  c  I  | 

13  kxk 


k 

where  [n]  **  I  n  ,  and  1 1  c  .11  is  the  kxk  determinant  with 
i«l  1  kxk 

(i,j)th  element 


_ 1 _ 

Cij  "  -  a  -  i  +  j)'.  * 

Kreweras'  original  result  is  in  terms  of  Young  chains  which 
when  converted  into  lattice  paths  has  the  above  interpretation. 

Before  stating  the  next  proposition  we  need  some  definitions 
and  notations.  Consider  (k  +  1) -dimensional  lattice  paths  from  the 
origin  to  the  point  (nQ,  n^...,^).  By  the  j:(=(r^, . . .  ,rfc))  th  level 

we  mean  the  set  of  points  {(x  ,  n,  -  r, , . . .  ,n.  -  r,  )  :  0  <  x  <  n  }. 

01  1  k  k  —  0—0 


x(r)  =  n  -  min{x  }  if  the  path  reaches  the  rth  level 

—  o  o 

and  (xo,  nx  -  ^ . -  rfc)  is 

a  point  on  the  path; 

-  0  otherwise . 

Clearly  0  <  x(r)  <  nQ  and  x(n)  **  nQ.  Denote  by  a(r)  and  b(r) 

the  upper  and  lower  restrictions  at  the  j:th  level,  by  which  we  mean  the 
path  at  the  rth  level  can  pass  only  through  points  in  the  set 

{(xq,  nx  -  ^,.-.,1^  -  rfc):  0  <  b(r)  <  nQ  -  xq  <  a(r)  <  nQ}. 

The  sets 

A(n)  *  {a(r)  :  ()  <  r^  n } 

and 

B(n)  -  (b(r):  £  <  v  <  n} 

are  respectively  called  the  upper  and  lower  restrictions  on  the  path.  The 
order  relation  x  <_  means  x^  <_  yi  for  each  i.  Note  that  a(r)  and 

b(r)  are  non-negative  integers,  and  non-decreasing  in  each  coordinate. 

Let  jE  a  2.  1116311  <  y±  for  at  least  one  i.  For  example, 

the  lexicographic  ordering  (u^, .. .  ,Uj)  of  the  set  {_r:  0.  <  r.  <  a} 
k 

such  that  d  *  TT  (n.  +  1),  u,  ■  _0  and  u,  •»  n  is  an  a  ordering  in 
i-1  1  1 

the  sense  that  u^a. ..au^.  Remember  that  the  sequence 

((0,0),  (0,1),  (1,0),  (1,1),  (2,0),  (2,1))  is  a  lexicographic  ordering 

of  vectors  (r:  (0,0)  <  <  (2,1)}. 

Proposition  2.  (Handa  and  Mohanty  (1979)) 

Denote  by  g^(A(n) lB(n))  the  number  of  paths  with  upper  re¬ 
striction  A(n)  and  lower  restriction  B(n) .  Let 


-  5  - 


k 

a(a  -  l)...(a  -  Z  n  +  1) 

_ i-l 

n! 


and 


k 

nl  *  TT  n  ! . 
i-1  1 

Then  g.  (A(n) | B ( n) )  satisfies  the  recurrence  relation 

(2)  2:  (-1)  “  I  r(-  +  V  g,  (A(r)|B(r»  -  6  J 

J)<r<n  V.  —  —  /  — 

where 

6^-1  when  n  «  0  , 

o  - 

-  0  otherwise. 

An  explicit  solution  of  (2)  is  the  following: 


d-[n]-l 


(3)  gk(A(n)lB(n))  »  (-1) 

k 


i>  -  t<iij+i>  A 
V  Sj+l  "  %  /+ 


where  d  - 


(d  -  1)  x  (d  -  1) 

17  (n  +  1)  and  {u. . u,}  »  {r:  0  <  r  <  n}  such  that 

i-1  1  _1  - 


0  -  u,  a . .  »au ,  -  n . 

—  — 1  — u  — 

The  explicit  expression  (3)  may  be  obtained  by  first  using 
Cramer's  rule  to  the  system  (2)  of  linear  equations  and  then  simplifying  it 


For  practical  purpose,  the  lexicographic  ordering  of  vectors  is  good 
enough.  It  is  easily  seen  that  if  we  take  (nQ,  n^,...,^)  as  the  origin 


and  reverse  the  steps  in  the  path,  the  upper  and  lower  restrictions 
respectively  become 


A'(n)  -  {nQ  -  b(n  -  jr) :  £  <  x  <  n } 

and 

B*  (n)  *  {nQ  -  a(ja  -  _r )  :  0  £  _r  <  n}  . 

Then  we  get  an  alternative  expression  for  the  same  number  of  paths  as 
*k(A’  (n)|B'(n)).  In  our  applications,  the  alternative  expression  is 
used. 

Next  we  formulate  a  continuous  analogue  of  Proposition  2  which 
gives  a  generalization  of  Steck's  result  (1971).  Though  Steck  has 
stated  the  result  in  terms  of  order  statistics,  its  connection  with 
paths  is  explained  In  Mohanty  (1979)  chapters  2  and  4;  Mohanty  (1980) . 

In  (k+1) -dimension  with  axes  xq,  x^,...x^,  consider  paths  (not 

necessarily  lattice  paths)  from  the  origin  to  (nQ,  n^,...,n^)  where 


nQ  is  a  non-negative  real  number  and  n^, . . . ,n^  are  non-negative  in¬ 
tegers.  In  this  case,  a  path  is  like  a  lattice  path  except  that  the 
the  number  of  units  moved  at  any  time  on  XQ-axis  is  a  non-negative  real 
number.  As  before,  we  may  define  the  level  £,  the  upper  restriction 
A(n)  and  the  lower  restriction  B(n).  Here,  we  may  remember  that 
a(r)  and  b(r)  are  non-negative  real  numbers.  In  the  next  assertion, 
we  adopt  some  of  the  earlier  notations. 

Proposition  3. 

Let  g^ACn)  |B(ja))  be  the  measure  of  the  set  of  paths  with 
upper  restriction  A(n)  and  lower  restriction  B(n)  .  Then 
g£(A(n)|B  (n))  satisfies  the  recurrence  relation 


a V ^  .It  !■!»  -Ii  -'a  Jj 


■  m.M  m  »  ■ 


\\\\ :  ii  i  uj l 


(4) 


(n  -  r)  i 


I  (-1) 

0  <  r  <  n 

max(0 ,x) 


gkU(£)  |B(r)) 


tn 


where  (x)  + 


An  explicit  solution  of  (4)  is  given  by 


(5)  g£(A(n)|B(n)) 


(-i)4"1— 1"1 


<Vi  -  Si>  ’• 


(d  -  1)  x  (d  -  1) 


The  proof  is  inductive,  follows  the  similar  line  as  in 
Proposition  2  and  therefore  is  omitted.  When  k  =  1,  we  get 


(6)  I  (-1) 
r=o 


n-r 


(a(r)  -  b(n)) 


n-r 


(n  -  r) 


y-i  gl(A(r)|B(r)>  -  6® 


and 


(7)  g. 


Cuj+1  '  ui^ 


(aC^)  -  b(Uj+j))+ 


<uj+i  *  V  '• 


n  x  n 


(g;L  written  for  g;L(A(n)  |B(n)))  . 

* 

In  fact  u^  *  r  -  1,  r  =  l,...,n  +  1  and  g^  represents  the  integral 


where  y^  =  max(u^,  x^  (see  Mohanty  (1971))  such  that 


a(r  -  1)  “  vr  a11'!  b(r)  =  u^  r  =*  l,...,n.  Thus  the  determinant  sol- 
* 

ution  (7)  of  g^  checks  with  the  earlier  one  in  Steck  (1971)  (see 
Mohanty  (1979),  page  56,  Mohanty  (1980)).  Relation  (6)  which  is  needed 
in  the  inductive  proof  follows  easily  from  (7). 
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The  remarks  following  Proposition  2  are  all  valid  for 
Proposition  3.  For  computational  purpose,  recurrence  relations  (2) 
and  (4)  are  often  more  useful  than  the  explicit  expressions  (3)  and  (5). 

3.  A  Tandem  Queue 

The  queueing  model  called  Model  I  is  described  by  the 
following  properties: 

Model  I. 

(a)  There  are  r  service  counters  numbered  1  to  r  and  each  counter  has 
only  one  server.  A  customer  arrives  at  the  first  service  counter  and 
moves  from  the  ith  counter  to  the  (i  +  l)st  counter  i  =  l,...,r  -  1 
for  service.  The  service  is  completed  just  after  the  customer  leaves 
the  rth  counter. 

(b)  Customers  are  arriving  in  accordance  to  a  Poisson  process  with 
parameter  X. 

(c)  Service  times  at  the  ith  counter  are  i.i.d.  exponential  random 
variables  with  parameter  y  ,  i  ■  l,...,r. 

(d)  Service  times  at  various  counters  are  mutually  independent  and  are 
independent  of  arrival  times. 

The  busy  period  ends  when  for  the  first  time  the  first  counter 
is  empty  while  the  other  counters  are  busy  and  are  never  empty  before. 

We  are  interested  in  ascertaining 

G^  -  G^(j^ . jf;  k^,...,kr_^;  n;  t)  which  represents  the  probability 

that  (i)  there  are  initially  j  +  1  customers  at  the  ith  counter, 
i  ■  l,...,r;  (ii)  n  customers  have  arrived  during  the  busy  period; 

(iii)  at  the  end  of  the  busy  period  there  are  k^  +  1  customers 
including  one  being  served,  left  at  the  ith  counter,  i  =  2,...,r; 


(iv)  the  length  of  the  busy  period  is  <t. 

Let  the  point  of  time  when  either  there  is  an  arrival  or  a 
completion  of  a  service  be  called  an  instant.  Denote  by  and 

0^(1  *  l,...,r)  the  number  of  customers  arrived  at  instant  j  and  the 
number  of  units  served  by  the  ith  server  at  instant  j  respectively.  Ob¬ 
serve  that  <3^  “0  or  1,  i*0,  l,...,r  such  that 


r 

Z  a. .  «  1.  Let  N  +  1  be  the  total  number  of  instants  upto  the  end  of 
i-o  1J 


the  busy  period.  Then  the  event  in  question  with  conditions  (i) ,  (ii) 
and  (iii)  is  equivalent  to  the  following  relations: 


N+l 


£-1 


£=1 


U 


ai£’  1  "  m  ■  1»***>N» 


l.N+l 


1. 


n  +  I  (j  -  k  )  +  k  +  1,  i  -  2,...,r. 
t-1 


ji+l+“,+jr 


m 

I 

£■ 


i,£  *i,m’ 


0,1,. 


r  -  1,  m  *  1 . N, 


m 

Z  a 
£-1 


ft  A 

r,£  r,m 


and 


(x  ,  x,  ,...,x  )  -  (I  j  ,  I  i . i  ,0) . 

oo  lo  ro  t~l  t  t=2  C  r 


we  note  that  relations  (9)  become 


10 


where 


%  =  ^XoN’  XlN’“*’XrN^ 


with 


VJ  ■  *  +  X1N  '  +  “■  “d  XiN 


Xo!l  '  ^kt  +  l’ 


X  3  2| 

The  number  of  sequences  ((x  ,  x.  _ ),  m  =  0,1,..., N) 

om  in  m 

satisfying  the  condition  xom  1  xlm  >  •••  >  xm  m  =  0,1,..., N  is  the 

number  N(r  +  1;  x^,  x^)  of  paths  as  in  Proposition  1. 

Taking  the  superposition  of  r  +  1  Poisson  processes  into 

consideration,  one  observes  that  there  are 

*  r  i  r 

n  *  (r  +  l)n  +  j  +  I  I  (j  -  k  )  +  I  k  +  r  -  1  occurrences  in  the 

1  i-2  t=l  C  i=2  1 


combined  Poisson  process  with  parameter  X  +  Z  p.  during  the  busv 

J-l  J 

period  of  length  y  and  one  occurrence  at  the  end  of  the  period  (i.e.,  dur¬ 
ing  the  interval  (y,  y  +  dy) ,  dy  ■+■  0)  .  The  probability  of  this  event  is 


(10)  e 


-y(X  +  Z  p)  (y(X  +  Z  p.))n 
1*1  J  -»=i  J 


((X  +  Z  p.)dy  +  ody) . 

j=l  J 


Given  the  above  event,  the  probability  for  any  sequence  of  arrivals  and 

completion  of  services  at  all  counters  during  the  busy  period  (i.e.,  any 

sequence  ((x  ,  x^  ,...,x  ),  m  =  0,1, ...,N  subject  to  the  condition 

om  1m  rm 


x  >  x.  >  ...  >  x  ,  m  =  0,1, .. . ,N)  is 
om  —  Xm  —  —  rm 


where 


o  r 

A  +  I  y , 


*  Pi  = 


3-1  J 


r 

A  +  I  y 


>  1  * 


3-1 


3 


Thus  coobiniag  these  facts,  we  obtain 


t  ~y(A+Zy.) 

C12)  Gx  -  N(r  +  li  r,,  %)  A  !o  .  J 


r 

where  £  refers  to  I  • 

1 


* 

(y.a  +  ry.))n 

' - — Z -  (A  +  Eyjdy 

n*» 


4 .  Markovian  queues 


We  deal  with  two  heterogeneous  Markovian  models  involving 
batches.  With  the  aid  of  Proposition  2,  the  joint  distribution  of 
several  characteristics  one  of  which  is  the  maximum  queue  length  during 
the  busy  period  is  obtained  for  each  model. 

Model  II  . 

The  model  is  characterized  by  the  following: 

(a)  Customers  arrive  at  a  counter  from  r  sources  in  batches  of  size 
from  the  ith  source  (i  =  l,...,r)  in  accordance  to  the  Poisson  process 

with  parameter 

(b)  Customers  are  served  in  batches  of  size  u  and  service  times  at 

o 

the  counter  are  i.i.d.  exponential  random  variables  with  parameter  y . 

(c)  Service  times  are  independent  of  arrival  times  from  different 
sources.  Arrival  times  from  different  sources  are  also  independent. 

Without  loss  of  generality,  we  will  first  treat  the  case 
r  =  2.  Let  +  uq;  n^  n2;  k,  t)  be  the  probability  that  a 

busy  period  initiated  by  l  +  uq  customers,  consists  of  n^  batches 

arriving  from  the  ith  source  (i  =  1,  2)  has  a  maximum  queue  length  <  k 
and  has  length  ^t. 

Let  nQ  +  1  be  the  number  of  batches  served  during  the  busy 
period.  Then 


u 

o 

[z]  being  the  largest  integer  less  than  or  equal  to  z.  Assume  that 

k  >  u  . 

—  o 


(13) 


V] 


-**■*-•  -*«  -*v*l*»^*»  -*•  -*.  -*-  -**  .*»  .vCVis  .**-.**  'j  '-i  'IuS  '■»  .'.t'1.  •-»  ' 


Let  us  define  the  (i,j)th  arrival  instant  as  that  instant  which  is 
either  the  ith  arrival  instant  from  the  first  source  and  j  arrival 
instants  from  the  second  source  precede  it  or  the  jth  arrival  in¬ 
stant  from  the  second  source  and  i  arrival  instants  from  the  first 


source  precede  it.  Denote  by  X^j  the  number  of  departures  after  the 


(i,j)th  arrival  instant  and  before  the  next  arrival  instant.  Also 


use  X 


al’n2 


to  represent  the  number  of  departures  between  the 


(n^,  n^Jth  arrival  instant  and  the  end  of  the  busy  period.  Then  for  a 


busy  period  the  following  relations  must  be  satisfied: 


i  i 

(14)  £  -  k  +  iu  +  ju_  <  u  I  X  X  <  iu  +  ju0  +  £ 

1  1  ~  °  &=o  a=0  1  2 

for  i  *  0,  l,...,n^>  j  a  0,  1,...,^  but  (i,j)  ^  (n^,^) 


n2  “l 

and  X  X  X  _  *  n  +  1. 

B-0  o-O  °’B  ° 

If  we  represent  a  departure  by  a  unit  on  xQ-axis,  an  arrival  from 
source  I  on  x^-axis  (i  ■  1,  2),  then  the  sequences  of  arrivals  and 
departures  satisfying  (14)  correspond  to  the  set  of  lattice  paths  from 
(0,  0,  0)  to  (nQ  +  1,  n^,  which  do  not  cross  the  planes 

(15)  u  x  ■  u,  x.  +  u_  x„  +  £ 
o  o  11  2  2 


and 

(16)  u  x  ■  u.  X.  +  u-  x„  +  £  -  k 
O  O  11  2  2 

except  at  the  end  when  the  plane  (15)  is  crossed  by  a  unit  step  on 
x  -axis.  The  number  of  such  paths  as  given  in  Proposition  2  is 


g2(A(n) | B (  n) ) ,  where 


iu.  +  ju_  +  £ 

a(i,j)  =  min(ao,  [ - - - ]) 

o 


£  -  k  +  iu.  +  ju, 
r - ± - 

b(i,j)  -  max(0,  uq  ) , 

i  ■  0,  l,...,n,  and  j  “  0,  1 . n2» 

Using  an  argument  similar  to  the  last  part  of  the  previous 
section,  we  obtain 

n  +1  n..  n.  t 

(17)  G2  =  g2(A(n)lB(n))  pQ°  p^  p^  /  B  dy 

o 

for  k  >  uq  and  nQ  given  by  (13) 
where 


U  +  ^  +  X2’  P1  y  +  XJL  +  X2  »  p2  y  +  +  X, 


n  +n,+n„ 


B  =  e 


-(y+X-+X,)y  «lri*X1+X2)y)  o  “1  “2 

1  2  7 -  --  -r,  (y  +  X.  +X). 

(n  +  n,  +  nil  1  2 


o  1 


When  k  <  uq,  the  busy  period  terminates  with  one  service  and  the  ex¬ 
pression  is  obvious. 

For  the  general  case,  we  may  write 

G^1^  ■  G2(£  +  uq;  n^,...,nr;  k,  t)  to  represent  the  probability  that 

busy  period  initiated  by  £  +  uq  customers  consists  of  n^  batches 
arriving  from  the  ith  source,  has  a  maximum  queue  length  <k  and  has 
length  <t.  The  value  of  n  becomes 


15  - 


Zu.n.  +  Z 

n  =  [-i-i - ]  . 

o  u  . 


The  expression  for  the  generalized  G2  is 


n  +1  r  n  t 

(17)  G2(r)  -  gr(A(n)  !B(n))po°  (TT  pii)  J  B  dy 

i*l  o 

where 

Is  u  +  l 

a(s)  *  min(n  ,  [ — - ]) 

—  o  u 


Is  u.  +  Z  -  k 

b(s)  =*  max(o,  [ - - ]) 


u 


for  0.  <  .s  <  n , 


li 


po  "  y  +  ZA±  *  p j  *  y  +  Z\± 
and 


»  ^  “  1» .  •  •  »r» 


b  .  .-<•*«!*  “u 

r  (n  +  Zh.)  : 

o  i 


n  +Zn . 
o  i 


(y  +  IX±) . 


In  the  above  discussion,  Z  represents  Z  . 

1 


Model  III. 

We  consider  another  Markovian  model  which  is  described  below: 

(a)  Customers  arrive  at  a  queueing  system  in  batches  of  size  u  in 

o 

accordance  to  a  Poisson  process  with  parameter  A. 

(b)  The  service  system  consists  of  r  counters.  Customers  are  served 
in  batches  of  size  u^  at  the  ith  counter  (i  *  l,...,r).  The  service 
times  at  the  ith  counter  are  i.i.d.  exponential  random  variables  with 
parameter  y^(i  *  l,...,r). 

(c)  Service  times  at  various  counters  are  mutually  independent  and  are 


Independent  of  arrival  times. 


The  busy  period  ends  when  for  the  first  time  any  one  counter 
becomes  empty*  and  without  loss  of  generality  let  us  assume  it  to  be  the 
last  counter. 

Because  of  our  observation  in  Model  II,  we  will  only  deal  with 
the  case  r  *  2.  Denote  by  +  u2*  n0»  nl* 

probability  that  a  busy  period  initially  with  l  +  u^  +  u2  customers 

consists  of  n  arrivals  of  batches  of  size  u  and  nn  batches  served 
o  01 

at  the  first  counter,  has  a  maximum  queue  length  <k  and  has  length 


Letting  +  1  to  represent  the  number  of  services  at  the  sec¬ 
ond  counter  during  the  busy  period,  it  can  be  verified  that 
n  u  +  Z  -  n,un 

U8)  n2  -  [  -°-°  - - —  ]  . 

2 

For  n  ,  n,  and  n_  to  be  meaningful,  we  assume  k  >  maxCu, ,  u„)  . 

O  i  4  —  Lt. 

As  in  Model  II,  we  define  the  (i,j)th  departure  instant  as 
that  instant  which  is  either  the  ith  departure  instant  at  the  first 
counter  and  j  departure  instants  at  the  second  counter  precede  it  or  the 
jth  departure  instant  at  the  second  counter  and  i  departure  instants 
precede  it.  Without  ambiguity,  we  may  again  denote  by  X^  the  number 
of  arrivals  after  the  (i,j)th  departure  instant  and  before  the  next 
departure  instant.  The  busy  period  will  consist  of  nQ  arrivals,  n^ 
batches  served  at  the  first  counter  and  n2  +  1  batches  served  at  the 


second  counter  when 

j  i 

(19)  u.i  +  u  j-i<  u  Z  Z  X  <  U]i  +  u2j 
A  A  p=o  a=o 


-  i  +  k  for 


wSVi 


such  that 


n2  Q1 

Z  Z  X.  *  n  . 

&-o  a-o  ^  ° 

Representing  an  arrival  by  a  unit  on  x^-axis  and  a  departure 
from  the  ith  counter  by  a  unit  on  x^-axis,  i  ■  1,  2,  the  number  of  se- 
quences  of  arrivals  and  departures  satisfying  (19)  can  be  seen  to  be 
equal  to  g2(A*(n) jBX(n))  where 

*  u.i  +  u-j  -  l  +  k 

a  (i,j)  -  min(nQ„  [— - - - ]), 

o 


*  M  +  “  1 

b  (i,j)  -  max(0,  [ - - - ))  , 

o 

i  "  0,  1  ,  •  •  •  , n^  and  j  °  0?  1 , •  •  •  1  • 

Again,  following  the  rest  of  the  argument  routinely  we  derive 

*  *  n  n1  n?+l  * 

(20)  G3  -  g2(A  (n)  '|B  (n))PQ°  P2  J  B  dy  for  k>max(u1,  i^) 

j  ,e.  O 


and  n2  given  by  (18) 


where 


/V  _ A _  A 

A  +  yx  +  ^2*  pl  "  X  +  +  y2’  p2  “  X  +  +  y2 


*  -(X+y1+y2)y  ((X  +  y1  +  y^y) 

B  "  6  <nQ  +  n1  +  n2)  ! 


no+nl+n2 


(A  +  ^  +  ^) 


Note  that  when  k  <  max(u^,  u2) ,  several  possibilities  arise. 
If  k  <  min(u3#  u2),  one  gets  n^  ■  n2  »  0.  On  the  other  hand. 


18  - 


■  0  if  u2  >  ui  and  <  k  <  u^  and  n^  *»  0  if  >  Uj  and 

Uj  <  k  <  u^.  Discussion  for  the  first  counter  to  be  empty  is  similar. 

The  derivation  for  the  general  case  is  similar  and  is  omitted. 
It  may  be  noted  that  the  models  in  Mohanty  (1972)  are  special 
cases  of  the  models  discussed  in  this  section. 
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5 .  Non -Markovian  Queues 

Two  models  in  which  either  the  input  or  the  service-time 
distribution  is  non-Markovian  are  considered.  Setting  inequalities 
similar  to  the  last  section  in  terms  of  either  arrival  instant  or  de¬ 
parture  instant,  an  expression  of  the  joint  distribution  involving  the 
nm-irfmn  workload  in  the  first  model  and  the  maximum  queue  length  in  the 
second  is  derived  with  the  help  of  Proposition  3.  The  advantage  of  this 
approach  is  mostly  achieved  when  the  non-Markovian  distribution  is  de¬ 
terministic. 

Model  IV. 

The  model  is  the  same  as  Model  II  except  that  (b)  is  changed 
to  the  following: 

(bf)  Customers  are  served  individually  and  the  service  times  {ij;  } 
of  customers  from  the  ith  source  are  i.i.d.  random  variables  with  distri¬ 
bution  function  H^(t)  i  »  l,...,r.  Service  times  for  each  source  are 
independent . 

Consider  only  the  case  r  -  2.  Let  G^  =  .(^(y;  n^  n2;  z,t) 

be  the  probability  that  a  busy  period  with  initial  work  load  y  con¬ 
sists  of  n^  batches  arriving  from  the  ith  source  (i  *  1,  2) ,  has  a 
maximum  work  load  <z  and  has  length  <t. 

Denote  by  f,  .  the  (i,j)th  arrival  instant  (see  Model  II  for 

J 

the  definition)  and  by  T^  the  instant  of  the  next  immediate  arrival 
instant.  A  busy  period  with  the  given  parameters  in  must  satisfy 

the  following  inequalities: 


E 


But  for  all  (i,j)  .  is  related  to  z  as  follows: 


y  -  T  +  I  +  I  ij,  <  z. 

X1  a-1  ia  a«l 


Therefore  (21)  becomes 


0  <  T  <  y 

—  oo  — 


lul  JU2  lul  3U2 

7  *  1  +  L,  *la  +  L,  *2a  i  TiJ  i  y  +  %  +  1  ,  * 


o-l  “  a-l  2a 


i  ■  0, j  m  0 , . . . , n2 »  (i, j)  ¥  (0»0). 


(Notice  the  similarity  between  (14)  and  (22)  .) 

Under  the  condition  that  there  are  n^  arrivals  from  the 
ith  source  (i  *  1,  2)  during  a  busy  period  of  length  t,  T 's  form 
two  Independent  sets  of  order  statistics  from  the  uniform  distribution 
over  (0,t),  the  joint  p.d.f.  of  which  is 
n* !  n, 

(23)  '  n  4ri 

nl^2 


Represent  the  arrival  time  on  XQ-axis  followed  by  a  unit 
on  x^-axis  if  the  arrival  is  from  the  ith  source.  Then  the  sequences 
of  arrival  instants  of  a  busy  period  are  represented  by  a  set  of  paths 
in  Proposition  3  each  of  which  has  the  weight  (23) .  Given  the  length 
of  the  busy  period  to  be  x,  the  upper  and  lower  restrictions 
A^(n)  and  B^(n)  of  the  set  are  given  by 

iul  iu2 

a,(i.j)  -  min(x,  y  +  I  +  I  %  ) 

a-1  xa  a-1 


and 

1U1  iu2 

b.(i,j)  -  max(0,  y  -  z  +  L  ip  +  Z  ip  ) 

a-1  X  a-1 

i  -  0,1,  ...,n^  and  j  -  0,1, 

Therefore,  given  there  are  n^  batches  of  arrival  from  the  ith 
source  (i  •  1,  2)  and 


“l  u2 

(24)  y  +  Z  *la  +  I  *2a  - 

a-1  a-1 

the  probability  of  the  busy  period  is 


(25) 


n  !  n  1 

X 


where  stands  for  the  conditional  expectation  taken  over 


and  subJect  to  the  condition  (24) 


Now  it  is  a  matter  of  routine  argument  as  before  to  establish 


22 


“l  P  2 

t  -(A  +A.)x  n  +n  A-  A,  ,.v  (2\ 

(29  0  -/  Ce  1  2  x1  2  i.A-  d(H<>(x-y)«H<>(x-y 

y  12  11  ^  ^ 

where  is  the  nth  convolution  of  with  itself  and  *  repre- 

n 

sents  the  convolution  operator. 

f  •  \ 

For  any  general  service  time  distribution  H''  '  expression 

(26)  seems  to  be  unwieldy.  However,  if  the  service  times  are 

deterministic,  namely,  the  service  time  for  the  customer  from  the  ith 

source  is  a  (i  *  1,  2) ,  then  the  probability  that  a  busy  period 

having  the  same  specifications  as  in  except  that  it  has  length 

t  -  y  +  a^u^n^  +  £-2u'>n2  is 

*  -(A.+A.,  )t  n  n 

(21)  g2(A1(n)l31(n))  e  1  1  \  X2 


where  the  specialized  values  of  restrictions  are  given  by 
a-jU.j)  »  min(t,  y  +  a^i  +  a^j) 


and 

b^Ci.j)  *  max(0,  y  -  z  +  a^i  +  a^j)  . 

Model  V. 

It  is  characterized  by  Model  III  except  that  (a)  is 
modified  as  follows: 

(a1)  Customers  arrive  one  at  a  time  at  a  queueing  system.  Let  <J>n 

denote  the  arrival  instant  of  the  nth  customer.  The  interarrival  times 

{d>  -  $  ,}  n  *  1,  2,...  are  i.i.d.  random  variables  with 

n  n-1 

distribution  function  M(t)  . 


Assume  that  the  busy  period  ends  when  the  first  counter  is 


empty.  For  r  -  2,  let  G5  -  G5(£  +  ^  +  u^.^  +  1,  n2;  k,  t)  be 

the  probability  that  a  busy  period  initially  with  Z^  +  ^  +  u2 

customers  consists  of  serving  n  +  1  batches  <.t  the  first  counter, 

Q 2  batches  at  the  second,  has  a  maximum  queue  length  <k  and  has 
length  <t . 

Without  any  confusion,  we  may  denote  by  x^  the  (i,j)th 
departure  instant  (see  Model  in  for  definition)  and  by  T  the  in- 

ij 

stant  of  the  next  immediate  departure  after  the  (i,j)th  departure  in¬ 
stant.  A  busy  period  with  the  constraints  of  G5  satisfies  the 
following  inequalities  t 

“0  <  T  <  d,  ,,  . 

—  oo-'Tc+l-f, 


(28) 


Tij-  Tijl  *iuj+  jUj-Ht+1-4 

|_fere 


for  i  -  0,  1,...,^,  j  -  0,  l,...,^,  (i,j)  $  {(0,0),  (nv  nj}  . 

If  we  fix  the  length  of  the  busy  period  to  be  equal  to  x,  we  have 


n.j+1,  n2 


V  n2 


x. 


Under  the  condition  that  there  are  n^  +  1  departures  from 
the  first  counter,  n2  from  the  second  during  the  busy  period  of  length 
x,  the  joint  p.d.f.  of  T^’s  is 


Represent  the  departure  time  on  x^-axis  followed  by  a  unit 
on  x.-axis  if  the  departure  is  from  the  ith  counter.  The  sequences  of 
departure  instants  become  a  set  of  paths  in  Proposition  3,  the  upper 
and  lower  restrictions  A^n)  and  B2(n)  being 


and 

b2(i,j)  =  max(0,  ^iUi+ju2_£) 
i  *  0,  l,...,n^  and  j  =  0,  l,...,n2> 

Therefore,  given  there  are  +  1  batches  of  departure  from  the 

first  counter  and  batches  from  the  second  and  the  busy  period 

is  of  length  x,  the  conditional  probability  of  the  busy  period  is 

n  T  „  t 

(29)  E*<82[<VH>IB2(2» 

X 

the  expectation  being  over  {$  }  subject  to  the  condition  that  the 
busy  period  is  of  length  x. 

The  busy  period  is  of  length  x  is  equivalent  to 
^  O  *  x,  where  K  =  n^u^  +  (n  2~  l)u  -  l  and  a  is  the  time  between 
the  Nth  arrival  and  the  last  departure.  Thus 


(30) 


n.j+1  n 
n,+n,  M,  h,/ 


G5  -  /  D  e"f  V^)X  x  1  2  n-1',-n--;  dt/X(l  -  M(x  -  u))d  J^(u)) 


1‘  2 


being  the  Nth  convolution  of  M  with  itself. 


When  the  interarrival  times  are  deterministic,  i.e. 


<j^  -  <+>n  ^  =  1  (say)  ,  the  probability  that  a  busy  period  having  the  same 
specifications  as  in  except  that  it  has  length  t  is 


-(y-j+y^t  n;L+l  n2 


yi  ^2 


(31)  g2  (A^lB^n))  e 

where  the  specialized  values  of  the  restrictions  are  given  by 
a2(i,j)  =  min(t,  iu_,  +  ju2  +  k  +  1  -  A) 

and 

b2(i,j)  =  max(0,  iu1  +  ju2  -  A). 
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